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Abstract 

This note extends the work of Foss and Tweedie (1998)| who showed 



that availability of the classic Propp and Wilson (1996) Coupling from The 



Past algorithm is essentially equivalent to uniform ergodicity for a Markov 
chain (see also Hober Tand Robert 20041) . In this note we show that all ge- 
ometrically ergodic chains possess dominated Coupling from The Past al- 
gorithms (not necessarily practical!) which are rather closely connected to 
Foster-Lyapunov criteria. 

1 Introduction 

Throughout this paper X will denote an aperiodic Harris-recurrent Markov chain 
on a measurable state space X which is a Polish space (the Polish condition is 
required in order to ensure existence of regular conditional probabilities). Recall 
that X is said to be geometrically ergodic if it converges in total variation and at 
geometric rate to statistical equilibrium n, with multiplicative constant depending 
on the starting point: 

dist TV (£(X n ),7r) < V(X )j n (1) 

for some function V : X — > [1, oo) and some rate 7 6 (0, 1). The chain X is said 
to be uniformly ergodic if the function V can be chosen to be constant. 
We also recall the notion of a small set: 
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Definition 1 A subset C C X is a small set (of order k) for the Markov chain X 
if there is a minorization condition: for (3 G (0, 1), and probability measure v, 

P [X k G E | X = x] > (3 1 [x G C] x z/(£) /or a// measurable EQX . 

(2) 

Results are often stated in terms of the more general notion of petite sets; however 
for ^-irreducible aperiodic chains the two notions are equivalent 
(Mey n and Tweedie 1993[ Theorem 5.5.7). 



poss and Tweedie (1998) use small set theory to show that the condition of 



uniform ergodicity for such X is equivalent to the existence of a Coupling from 
the Past algorithm in the sense of |Propp and Wilson (1996)} This classic CFTP al- 
gorithm delivers a perfect sample from the equilibrium distribution of X. The key 
to the F oss and Tweediel argument is to remark that in case of uniform ergodicity 
the entire state space is small. Sub-sampling the process X if necessary (to reduce 
the lorder of the small setl to 1), one can then devise a classic CFTP algorithm 
which is actually of the form introduced by |Murdoch~a nd Green (1998) as the 



multigamma coupler. Hobert and Robert (2004)| develop the IFoss and Twe edie 
argument to produce approximations to deal with burn-in (time till approximate 
equilibrium) in the geometrically ergodic case. 

The Foss and Tweediel result might be thought to delimit and constrain the 
possible range of applicability of CFTP. However it is also possible to sam- 
ple perfectly from the equilibrium of some strictly geometrically ergodic chains 
using a generalization: namely dominated CFTP (domCFTP) as introduced in 
Kendall (19981 IKendall and M0ller (2000)| |Cai and Kendall (2002)| In this note 



we show that this is generic: geometric ergodicity implies the existence of a spe- 
cial form of domCFTP algorithm adapted to the geometric ergodicity in question. 
Recent expositions of quantitative convergence rate estimation depend heavily on 
small sets and their relatives (see for example iRosenthal 20 02 ). so this piece of 
CFTP theory connects to quantitative convergence theory in a rather satisfying 
way. 

To describe this special form of domCFTP, we must first introduce the notion 
of a Foster-Lyapunov condition. Geometric ergodicity for our X is equivalent to a 
geometric Foster-Lyapunov condition involving recurrence on small sets (this can 
be extracted from |Meyn" and Tweedi e 1993[ Theorem 16.0.1): 

E [A(X n+1 ) \X n = x] < aA(x) + b I [X n G C] , (3) 

for some a G (0, 1) and b > 0, some lsmall setl C. and a function A : X — > [1, oo) 



2 



which is bounded on C. Note that a + b > 1 is required, as is A\ C c > a -1 , since 
we impose A > 1. 

Now Condition (3)| implies that every sub-level set {x E X : A(x) < c} is 



small (as indeed do weaker conditions; Meyn and Tweedi e~1993[ Theorem 14.2.3). 



This is a key fact for our argument so we sketch a coupling proof. 

First note that without loss of generality we can employ sub-sampling to en- 
sure that the small set C in Condition (3)| is of lorder ll Super-martingale argu- 



ments show that we can choose n such that P [X hits C before n \ X = x) can 
be bounded away from zero uniformly in x for A(x) < c. Let the hitting prob- 
ability lower bound be po- We can use the Minorization Condition (2)| to realize 



X as a split-chain in the sense of [Nummeli n (197%)] , regenerating with probabil- 
ity (3 whenever X E C. Couple chains from different starting points according 
to the time when X first regenerates in C, yielding a family of realizations X x 
of the Markov chain, with Xq = x, such that with positive probability (3p all 
realizations {X x : A(x) < c] coalesce into a set of at most n trajectories by 
time n (divided according to the time of first regeneration). Now apply a renewal- 
theoretic argument to the subsequent regenerations of this finite set of trajectories, 
which are allowed to evolve independently, except that whenever two trajectories 
regenerate at the same time they are forced to coalesce. Straightforward analysis 
shows that we can choose m such that with positive probability p\ < (3p all tra- 
jectories starting from {x E X : A(x) < c] have coalesced to just one trajectory 
by time n + m. Hence {x E X : A{x) < c] is a small set of order n + m, 
with minorization probability p x . It is convenient to isolate the notion of a scale 
function such as A in |Equation (3) 



Definition 2 A (Foster- Lyapunov) scale function for a Markov chain state space 
X is a measurable function 

A : X -> [l,oo) 

such that sub-level sets {x E X : A(x) < A} are small for all A > 1. 

Now we can define the special form of domCFTP which we require, which is 
adapted to a specified Foster-Lyapunov scale function. 

Definition 3 Suppose that A is a scale function for an Harris-recurrent Markov 
chain X. We say the stationary ergodic random process Y on [1, oo) is a dominat- 
ing process for X based on the scale function A (with threshold h and coalescence 
probability e) if it is coupled co-adaptively to realizations of X x ~ t (the Markov 
chain X begun at x at time —t) as follows: 
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(a) for all x G X, n > 0, and —t < 0, almost surely 

A(*-i+„) < Y-t+n => HX-t+n+l) < ^-t+n+l! (4) 

moreover ifY n < h then the probability of coalescence is at least e, where 
coalescence means that the set 

{X*£* : such that - t<n and A(X^~ t ) <Y n ) 

is a singleton set; 

(c) and finally, P [Y n < h] must be positive. 

Suppose Y is a dominating process for X based on the scale A. The following 
domCFTP algorithm then yields a draw from the equilibrium distribution of X. 

Algorithm 4 

Simulate Y backwards in equilibrium till the most recent T < for which 

Y T < h; 

while coalescence does not occur at time T: 

extend Y backwards till the most recent S < T for which Y$ < h; 
set T <— S; 

simulate the coupled X forwards from time T + 1, starting with the unique 
state produced by the coalescence event at time T; 

return X as a perfect draw from equilibrium. 

Practical implementation considerations are: (1) can one draw from the equilib- 
rium of y? (2) can one simulate Y backwards in equilibrium? (3) can one couple 
the dominated target processes X x, ~ t with Y so as to ensure the possibility of re- 
generation? (4) can one determine when this regeneration has occurred? and, of 
course, (5) will the algorithm not run too slowly? 

The simplest kind of ordinary small-set CFTP, as in Murd ocE and Green (1 998) , 
is recovered from this Algorithm by taking Y = h, and requiring the whole 
state-space to be small. In actual constructions, care must be taken to ensure 
that Y dominates a coupled collection of X for which coalescence is possible 
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as specified in [Definition 3(b)| (see the treatment of CFTP for Harris chains in 



ICorcoran and Tweedie 2001 1) . 

The proof that this algorithm returns a perfect draw from the equilibrium dis- 
tribution of X is an easy variation on the usual domCFTP argument, found at vary- 
ing levels of generality in Kendall 1991 |Kendall and M0ller ICai and Kendall 2(W\. 
The key is to observe that [Algorithm 4 reconstructs a coalesced trajectory which 
may be viewed as produced by the Markov chain begun at time — oo at some 
specified state x such that A(x) < h: the proof is then an exercise in making this 
heuristic precise. 

The poss and Tweedie (1998)| argument, and the fact that the geometric Foster- 
Lyapunov [condition (3)| would certainly produce a dominating process if the ex- 
pectation inequality was replaced by a stochastic domination, suggests our main 
result, which will be proved in lSection~2l 

Theorem 5 If X is a geometrically ergodic Markov chain, and A is a scale func- 
tion for X which is derived from some geometric Foster-Lyapunov condition, then 
there exists a domCFTP algorithm for X (possible subject to sub-sampling) using 
a dominating process based on the scale A, as in Algorithm~4\ 



As in the case of the |Foss and Tweedie ( 1998) result, this algorithm need not 
be at all practical! 

2 Geometric ergodicity implies domCFTP 

We begin with a lemma concerning the effect of sub-sampling on the geometric 
Foster-Lyapunov I condition! 



Lemma 6 Suppose X satisfies a geometric Foster-Lyapunov condition] for some 



a < 1, some scale function A, and small set C = {x G X : A(x) < c}. 

E[A(X n+1 ) \X n = x] < aA(x) + bI[A(X n )<c] . (5) 
Under k-sub-sampling we obtain a similar condition but with different constants: 

E [A(X n+k ) \X n = x] < a k - l A(x) + V I [A{X n ) < c'} , (6) 
and also, ifk>2, 

E [A(X n+k ) \X n = x] < aA{x) + b" I [A(X n ) < c"\ . (7) 
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Moreover b' = 6/(1 — a), cf = 6/ (a 1 (1 — a) 2 ) may be chosen not to depend on 
c, and b" = 6/(1 — a), c" = 6/ (a(l — a) 2 ) may be chosen to depend neither on c 
nor on k > 2. 

We are able to choose 6', el, b", c" not to depend on c because we have allowed 
generous sub-sampling (i.e.: /c-sub-sampling to change a to c^" 1 ). 
Proof: Iterating [Equation (5)[ 

k 

E [A(X n+k ) \X n = x] < a k A(x) + ^ a^bE [I [A(X n+k ^) < c] \ X n = x] 

3=1 

< a k A(x) H — 

1 — a 

„ii-i * /„\ „.k-ir- ■ b 



a«- l A(x) - a* _i (l - a)A(x) + 



1 — a 



Q ! fc - 1 A(x) if A(a?) > 



) a fc_1 A(x) + 6/(1 — a) otherwise. 

Hence we may choose 6' = 6/(1 — a), c' = 6/ (a fc_1 (l — a) 2 ). Alternatively 

6 

E \A(X n+k ) \X n = x] < aA(x) - ail - a k - l )A(x) + 

1 — a 

a H%) lf A (^) > a (i- a) (i- a *-i) ' 

aA(x) + 6/(1 — a) otherwise. 



< 



Hence we may choose b" = 6/(1 — a), c" = 6/ (a(l — a) 2 ) if k > 2. □ 

Proof (of Theorem 5): 

We first construct the dominating process. 



Consider Markov's inequality applied to the geometric Foster-Lyapunov inequality (3) 
Any dominating process Y must satisfy the |stochastic domination (4)| described in 
Definition 31 Consequently, in default of further distributional information about 
P [A(X n+ i) \X n = x), if Y is to be a dominating process based on the scale A then 
we need Y to be stationary ergodic but also to satisfy 

ply , lv ! . E[A(X n+1 )\X n = x] 

F [Yn+i > azy \Y n = z\ > sup . (8) 

x:A(x)<z OiZy 
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Now if C C {x e X : A(x) < c} then 

E [A(X n+1 ) I X n = x] aA(x) + bl [x : A(x) < c] 
sup < sup 

x:A(x)<z x:A(x)<z OCZJJ 

aA(x) 1 6 

< sup = - so long as z > c H — . 

x-.A(x)<z azy y " a 

Consequently Y is a possible candidate for a dominating process based on the 
scale A if 

ppu>™»in = *] = (; /!/ tf ** c+i ' c) 

I 1 otherwise. 

If we define U by Y = (c + b/a) exp(C7) (so C/ is a log-dominating process) then 
£7 is the system workload of a D/M/l queue, sampled at arrivals, with arrivals 
every log(l/a) units of time, and service times being independent and of unit 
Exponential distribution. The process U is a random walk with reflection (of 
Skorokhod type) at 0: as its jump distribution is Exponential 1) — log(l/a) we 
may deduce it is positive-recurrent if and only if a < e -1 . 

In case e -1 < a < 1, U and Y = (c + b/a) exp(U) fail to be positive- 
recurrent. However the same construction will work if we use Equation (6)| of 



ILemma 6l to justify sub-sampling X with a sampling period k large enough to 
ensure a geometric Foster-Lyapunov condition (3) using A as scale but with a re- 
placed by ct fc_1 < e _1 , and amending b to b', c to d as in [Inequality (6) 



Thus without loss of generality we may assume a < e , and so this Y can 
be run in statistical equilibrium, and thus qualifies as least partly as a dominating 
process for the purposes of ITheorem 51 In the sequel we assume moreover that 
further sub-sampling has been carried out based on Equation (7)[ to ensure that 
the following small set is of order 1: 

{x E X : A fx) < h} for h = max { c + -, — — ( 1 + ' 



a a(l — a) \ 1 — a 

(10) 

Here the level h > c + b/a is fixed so as to ensure h = c" + b"/(l — a) with b", c" 
given as in Equation (7)[ thus h supplies a stable threshold for geometric Foster- 



Lyapunov conditions, even allowing for further sub-sampling if required. Note in 
particular that Y — (c + b/a) exp(C7) is able to sink below h, since h > c + b/a 
and the system workload U can reach zero. 
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To fulfil the requirements on a dominating process given in Definition 31 we 
need to construct a coupling between Y and the target process X expressed in 
terms of a random flow of independent maps F- t +n+i '■ X -> X: 

v x,—t _ t-j i y x '~ * "\ 

satisfying the distributional requirement that X x ~ l should evolve as the Markov 
chain X, the domination requirement expressed by the implication (4)[ and also 



the regeneration requirement that with probability e the set 

{F n {u) : such that A (u) < h} 

should be a singleton set. The well-known link between stochastic domination 
and coupling can be applied together with the arguments preceding |Equation (9) 



to show that we can couple the various X x ~ l with Y co-adaptively in this manner 
so that the implication (4) holds: note that here and here alone we use the Polish 
space nature of X, which allows us to complete the couplings by constructing 
regular conditional probability distributions for the various X x ~ t conditioned on 
the A(X x ~ t ). Thus all that is required is to show that this stochastic domination 
coupling can be modified to allow for regeneration. 

The small set condition for {x E X : A(x) < h} means there is a probabil- 
ity measure v and a scalar j3 G (0, 1) such that for all Borel sets B C [1, oo), 
whenever A(x) < h, 

F[A(X n+1 )eB\X n = x] > Pv(B). (11) 

Moreover the stochastic domination which has been arranged in the course of 
defining Y means that for all real u, whenever A(x) < y, 

P [A(X n+1 ) > u | X n = x] < P [Y > u | Y = y] . (12) 

We can couple in order to arrange for regeneration if we can identify a probability 
measure v, defined solely in terms of v and the dominating jump distribution 
P [Y > u | Y = y], such that for all real u 

P [A(X n+1 ) >u\X n = x]- /3i/((it, oo)) < ¥[Y>u\Y — y]— {3v((u, oo)) 

v((u,oo)) < v((u,oo)) 

and moreover 
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For then at each step we may determine whether or not regeneration has occurred 
(with probability /3); under regeneration we use stochastic domination to couple v 
to v; otherwise we use stochastic domination to couple the residuals. 

We state and prove this as an interior lemma, as it may be of wider interest. 

Lemma 7 Suppose U, V are two random variables defined on [1, oo) such that 

(a) The distribution C (U) is stochastically dominated by the distribution C (V): 

¥[U>u] < F[V>u] for all real U ; (13) 

(b) U satisfies a minorization condition: for some (3 G (0, 1) and probability 
measure u: B C [1, oo), 

P [U E B] > /3u(B) for all Borel sets B C [1, oo) . (14) 

Then there is a probability measure \x stochastically dominating v and such that 
P/j, is minorized by C (V). Moreover \i depends only on (3v and C (V). 

Proof (of Lemma 0: 

Subtract the measure (3v{{u, oo)) from both sides of [Inequality (13)| representing 
the stochastic domination C (U) ^ C (V). By the minorization condition (14) the 



resulting left-hand-side is nonnegtive. Thus for all real u 

< P [U > u] - Pu{(u, oo)) < F[V >u\- (3v{{u,oo)) 



Now £ (U)—(3u is a nonnegative measure (because of the |minorization condition (14) I 



Consequently P [U > u]—(3v((u, oo)) must be non-increasing in u and so we may 
reduce the right-hand side by minimizing over w < u: 

¥[U > u] -/?!/((«, oo)) < inf {¥[V > w] - /3i/((w,oo))} 

w<u 

= ¥[V > u] -Pfi((u,oo)) 
where \x is the potentially signed measure defined by 

Pfi([l,u]) = W[V <u]-swp{F[V <w]- Pv([l,w))} . 

w<u 

In fact \i is a probability measure on [1, oo). Both/i({l}) = ^({1}) and/i([l, oo)) 
1 follow from considering u = 1, u — > oo. Now we show /x is nonnegative: 

Pn((u, u + u'}) - P [u < V < u + u'} 

= - sup {F[V < w] - (3v{[l,w))} + sup{F[V < w] - /3u([l,w))} . 

w<~u-\-u' W<U 



9 



If the first supremum were to be attained at w < u then the two suprema would 
cancel. If the first supremum were to be attained at w' E [u,u + u'} then 

Pn((u, u + u']) - P [u < V < u + u] 

= —¥[V<w'] + Pu([l,w')) + sup {P [V<w]- Pu([l,w))} 

w<u 

> -¥[V<w'] + pv([l, w')) + P [V < u] - pu([l, u) 
and hence 

Pn((u,u + u'}) > F[w' < V < u + u] +Pv([u,w')) > 0. 

So we can deduce (3fi is in fact a nonnegative measure. 
On the other hand 

Pfi((u, u + u']) - P [u < V < u + u'] 

= - sup {F[V < w] -Pv([l,w))} + sup{P[V < w] -Pv{[l,w))} 

w<.u-\-u' W<U 

< o, 

hence 

< Pfj,((u,u + u']) < F[u<V<u + u'), (15) 
so P/j, is absolutely continuous with respect to C (V) and indeed we can deduce 

Pdfji(u) = I [P [V > ■} - Pu((-, oo)) hits current minimum at u] dP [V < u] . 

(16) 

The minorization of Pfj, by £ (V) follows from this argument: dependence only 
on Pv and C (V) follows by construction; finally, stochastic domination of Pv 
follows from 

Pfi((u, oo)) = ¥[V > u] - inf {F[V > w] - pu((w,oo))} 

w<u 

= sup{/?z/((w,oo)) -F[w < V < u}} 

w<u 

> Pv((u, oo)). 

□ 

This concludes the proof of Theorem 5: use lLemma 7l to couple C (X n+1 | X n = x) 
to L (Y n+ i \Y n = y) whenever A(x) < y in a way which implements stochastic 
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domination and ensures all the X n+ \ regenerate simultaneously whenever Y < h. 

□ 

Note that the algorithm requires us to be able to draw from the equilibrium 
distribution of Y and to simulate its time-reversed equilibrium dual. Up to an ad- 
ditive constant log(Y) is the workload of a D/M/l queue. This queue is amenable 
to exact calculations, so these simulation tasks are easy to implement (specializ- 
ing the theory of the G/M/l queue as discussed in Gri mmett and Stirzaker 1 992, 
ch. 11). However in general we do not expect this "universal dominating process" 
to lead to practical domCFTP algorithms! The difficulty in application will arise 
in determining whether or not regeneration has occurred as in [Algorithm 4 This 



will be difficult especially if sub-sampling has been applied, since then one will 
need detailed knowledge of convolutions of the probability kernel for X (poten- 
tially a harder problem than sampling from equilibrium!). 

Of course, in practice one uses different dominating processes better adapted 
to the problem at hand. For example an M/D/l queue serves as a good log- 
dominating process for perpetuity-type problems and gives very rapid domCFTP 
algorithms indeed, especially when combined with other perfect simulation ideas 
such as multishift CFTP (IWilson 2(XX)b> . read-once CFTP (IWilson 20 00a). or 
one- shot coupling (Roberts and Rosenthal 20021 . 

Finally note that, in cases when a E [e -1 , 1) or when the small set {x E 
X : A(x) < h} is of order greater than 1, we are forced to work with coupling 
constructions that are effectively non-co-adapted (sub-sampling means that target 
transitions X mk to X mk+1 depend on sequences Y mk , Y mk+1 , Y mk+k ). The po- 
tential improvements gained by working with non-adapted couplings are already 
known not only to theory (the non-co-adapted filling couplings of Griffe ath" 19751 
Goldste in 19791 and the efficiency considerations of Burdzy and Kendall 2000) 
but also to practitioners (IHuber 20041 non-Markovian techniques in CFTP; 



|Hayes and Vigoda 2003 non-Markovian conventional MCMC for random sam- 



pling of colorings). 



3 Counter-example 

We complete this note by describing a counter-example: a Markov chain X which 
satisfies a Foster-Lyapunov condition involving a scale function A, but such that 
there can be no recurrent dominating process Y based on A. We begin by choosing 
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a sequence of disjoint measurable sets Si, S 2 , . . . , subsets of [1, oo) such that each 
set places positive measure in every non-empty open set: 

Lemma 8 One can construct a measurable partition Si, S 2 , ■ ■ . of [1, oo), 

SiUS 2 US 3 U... = [l,oo), 
with the property Leb(Si fl (u, v)) > for all < u < v < oo, all i G {1,2,.. .}. 

Proof: Enumerate the rational numbers in [0, 1) by = q , qi, q 2 , Choose 

a < 1/2, and define 

oo oo 

A o = \J{J[qn + k,q n + k + a2~ n ] . 

k=l n=0 

Then for each k > 1 

a < Leb(A n [k,k + 1)) < 2a. 
Continue by defining a sequence of nested subsets A r C A r _ x by 

oo oo 

Ar = UU 

k=l n=0 

satisfying 

a f . [ k k + 1\\ 2a 

— < LebUfl — , — — < — . (18) 

Thus the measurable shell B r = A r \ A r+1 places mass of at least in each 
interval [£, ^1) . 

It follows that if 5 is defined by 

oo 

S = \J(A rs \A rs+1 ) 

8 = 1 

then Leb(S' fl U) > for every open set [/ C [1, oo). The desired disjoint se- 
quence Si, S 2 , ... is obtained by considering a countably infinite family of disjoint 
increasing subsequences of the natural numbers. □ 



Qn + k q n + k 



a 



(17) 
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Lemma 9 There is a Markov chain X satisfying a Foster-Lyapunov condition 
with scale function A, such that any dominating process Y based on A will fail to 
be positive-recurrent. 

Proof: The Markov chain X will have state space [1, oo), with scale function 
A(x) = x. We begin by fixing a E (e _1 , 1), and set C = [1, a -1 ]. The set C will 
be the small set for the Foster-Lyapunov condition. Choose a measurable partition 
S% U S 2 U S3 U . . . = [1, 00) as in lLemma~51 Enumerate the rational numbers in 
[1, 00) by q x , q 2 , .... 

We define the transition kernel p(x, •) of X on [1, 00) as follows: 

For x E [1, a -1 ], set 

p{x,dy) = exp(-(y - l))dy fmy>l, 

so that if X n E C then X n+1 — 1 has a unit rate Exponential distribution. 
Then: 

C is a small set for X of order 1 (in fact it will be a regenerative atom!); 
if X n E C then E [X n+1 ] = 2; 

if X has positive chance of visiting state 1 then the whole state space 
[1, 00) will be maximally Leb-irreducible. 

For x > or 1 and x E Si, set 

(oc\ a 
1 5 (dy) + -5 mx (dy) . 

Note that, because we are using the identity scale A (a;) = x, 

ifx^C then E [A(X n+1 ) \ X n = x] = E [X n+1 \ X n = x] = ax; 
ifxE'C then P [X n+1 = 1 \ X n = x] > 0. 

Thus X satisfies a geometric Foster-Lyapunov condition based on scale A and 
small set C, and so is geometrically ergodic. 

Suppose Y is a dominating process for X based on the identity scale A. This 
means it must be possible to couple Y and X such that, if A(X n ) = X n < Y n 
then A(X n+ i) = X n+ i < Y n+ i. This can be achieved if and only if 

P [X n+1 >z\X n = u] < P [Y n+1 >z\Y n = x] 
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for all z > 1, and Lebesgue-almost all u < x. Therefore we require of such Y 
that 

P [Y n +i > axy \Y n = x] > ess sup {P [X n+X > axy \ X n = u]} 

u<x 

= sup ess sup < — : a^ 1 < u < x,u G Si, qiu > axy > 

i Ui J 

fa \ 1 

= sup <—: qi> ay > = -, 
i I ft J 2/ 

using Markov's inequality, then the construction of the kernel of X, then the 
measure-density of the Si. 

So such a Markov chain Y must also (at least when above level a -1 ) domi- 
nate exp(Z), where Z is a random walk with jump distribution Exponential 1) + 
log(a) . Hence it will fail to be positive-recurrent on the small set C when a > e _1 . 

□ 

There may exist some subtle re-ordering to provide domCFTP for such a chain 
on a different scale; however the above lemma shows that domCFTP must fail for 
dominating processes for X based on the scale A. 



4 Conclusion 

We have shown that geometric ergodicity (more strictly, a geometric Foster— 
Lyapunov condition) implies the existence of a special kind of domCFTP al- 
gorithm. The algorithm is not expected to be practical: however it connects 
perfect simulation firmly with more theoretical convergence results in the spirit 
of the F oss and T weedie '(1998)| equivalence between classic CFTP and uniform 
ergodicity. Note also that the "universal dominating process", the sub-critical 
exp(£)/M/l) so derived, is itself geometrically ergodic. 

It is natural to ask whether other kinds of ergodicity (for example, polynomial 
ergodicity) can also be related to perfect simulation in this way; this is now being 
pursued by Stephen Connor as part of his PhD research at Warwick. 
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